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Abstract. The form factor of the unitary group U(N) endowed with the Haar 
measure characterizes the correlations within the spectrum of a typical unitary matrix. 
It can be decomposed into a sum over pairs of "periodic orbits" , where by periodic orbit 
we understand any sequence of matrix indices. From here the diagonal approximation 
can be defined in the usual fashion as a sum only over pairs of identical orbits. We 
prove that as we take the dimension N to infinity, the diagonal approximation becomes 
"exact" , that is converges to the full form factor. 
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1. Introduction 

One of the central questions in the quantum theory of classically chaotic systems is 
whether the fluctuations in the spectra follow the predictions of random matrix theory 
(RMT). The spectral form factor, which is the Fourier transform of the two-point 
correlation function, is often used to characterize these fluctuations. 

The advantage of the form factor is that it has a convenient expansion in terms of 
pairs of periodic orbits of the classical system. This expansion is an application of trace 
formulae (see £Q for a review) and has been a starting point for many investigations. The 
progress of understanding the role of periodic orbits in the universal behavior of the form 
factor of different systems is marked by such milestones as the diagonal approximation 
[2], the first off-diagonal contribution j3Hj and the most recent preprint [5] outlining the 
complete expansion. Quantum graphs were also used [6 j as a mathematical toy model 
with a simple-to-understand periodic orbit structure and an exact trace formula. 

In a recent paper, Degli Esposti and Knauf \7\ proposed to use the ultimate 
mathematical toy model for the random matrix conjecture — the random matrix theory 
itself. The authors argued that if one can perform a step-by-step expansion of the RMT 
form factor in terms of families of periodic orbits, it would greatly facilitate analogous 
expansions for other systems. 

Degli Esposti and Knauf [Zj studied the diagonal approximation to the form factor 
of one the random matrix ensembles, the Circular Unitary Ensemble (CUE). Their 
results, however, depend on a combinatorial conjecture, which is still unproven. The 
purpose of this letter is to present an alternative proof for their main result; our proof 
is unconditional and is significantly simpler. Our proof is based on the observation that 
the diagonal approximation can be expressed as a trace of certain doubly stochastic 
matrix. The asymptotic properties of the diagonal approximation can then be related 
to the spectral gap of the doubly stochastic matrix, which was previously studied by 
the author in [T2] . 

The structure of the letter is as follows. In Section |2] we introduce the Circular 
Unitary Ensemble and the form factor. We follow with reviewing the main result of [7] 
and stating our main result in Section El The relevant properties of doubly stochastic 
matrices are reviewed in Section EJ Here we also prove our main result which follows 
quite simply from a technical lemma on the asymptotics of the second largest eigenvalue. 
This lemma is proved in Section El and we conclude with comparing our results with 
those of [7j in Section |H1 

2. Circular Unitary Ensemble and its form factor 

The Circular Unitary Ensemble is the unitary group U(iV) equipped with the Haar 
measure: 

Definition 1. CUE(iV) is defined as the ensemble of all unitary NxN matrices endowed 
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with the probability measure that is invariant under every automorphism 

U i-> VUW, (1) 

where V and W are any two N x N unitary matrices. 

The form factor is a function describing the statistical properties of the spectrum of 
U. It is defined as the Fourier transform of the two-point spectral correlation function, 



\ 3,k=l I N 

where exp(27rz0j) are the eigenvalues of U and the averaging is performed over the 
CUE(iV). The form- factor of the unitary ensemble is thus given by 

K N {t) = ±{\Tr(U t )\% (ieZ), (3) 
The well-known formula for K^{t) is 

N, t = 

K N (t) = { \t\/N, 0<\t\<N (teZ). (4) 
1, N<\t\ 

This formula can, in particular, be derived by expanding the trace in (J3J), 

TV TV / t \ 

ii,...,it=l ji,...,jt=l \k=l I TV 

where the convention i t+ i = %\ is assumed here and throughout the manuscript. The 
above expression is analogous to the trace-formula expansions of the form factor in other 
systems, with i = (ii, . . . ,i t ) and j = (ji, . . . ,j t ) playing the role of the periodic orbits. 
It can be shown (see equation JSJ) below) that only the terms with j being a permutation 
of i contribute to the form factor, i.e. 

Kn= Ii e e (n^iw^« + i)) ' ( 6 ) 

ii,...,tt=l aeSt \k=l I TV 

where St is the group of all t-permutations (the symmetric group). This step is analogous 
to only considering pairs of orbits with a small action difference. The approximation to 
(JHJ) that restricts the double sum even further, to pairs of periodic orbits that are equal 
(up to a cyclic shift), is called the diagonal approximation 

N I t 



k ^ = n E (lliiwr) > ( 7 ) 

il,...,it=l \k=l I TV 

where the factor t comes from the number of possible shifts of an orbit. In this letter 
we compare the form factor K^{t) with its diagonal approximation. 

The present article was inspired by a study by Degli Esposti and Knauf, [Zj. 
They proved that, assuming validity of a certain combinatorial conjecture, the diagonal 
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approximation converges to t/N uniformly in t/N as iV — > oo. The convergence 

is proved in a closed interval which is bounded away from zero. The combinatorial 
conjecture requires a more thorough introduction and will be delayed to a later section. 

Here we propose an alternative approach which exploits the fact that the sum in 
Ki t N is the trace of a doubly stochastic matrix. Our proof is unconditional, i.e. it does 
not rely on the conjecture, and the interval of convergence is extended. 

3. Averages over unitary group and the main result 

The averages of products of the elements of matrices U have been studied by various 
authors, [HI El UH1 [TTj . Their main result is 

t 

(u aibl ...u asbs u: iPi ...u: tPt ) N = 5 s t J2 M^iIIC'C". ( g ) 

cr,ir€St k=l 

where St is the symmetric group of permutations of the set {1, . . . , t}, 5% is the Kronecker 
delta and the coefficient Vn{c~ 1 'k) depends only on the lengths of cycles in the cycle 
expansion of o~ Xr n. 

To introduce the results of [7] we need the following definitions. 

Definition 2. Let a £ St and i £ {1,2, ... ,t}. Then a cycle of i with respect to o is 
the (finite) set {o~ k (i): k £ N}. The cycle length of i is the number of elements in the 
cycle of i. The rank \o~\ of a permutation a is the number of disjoint cycles induced by 
a. By \o~ V p\ we denote the number of disjoint sets of the form {p n a (i): k,n £ N}. By 
r we will denote the circular permutation (1,2, ... ,t). If a £ S t , o~+ will denote r^Vr. 

Conjecture 1 (Degli Esposti — Knauf |7J). There exists a constant C% such that for 
allt<NeN 

V N (tt~V) ATl^+l < C^iVl^+H a £ S t , (9) 

ireSt 

Theorem 1 (Degli Esposti — Knauf [7J). Assuming Conjecture^ is true, the form 
factor Kn is approximated by the diagonal contribution in the following sense: 
For all e > 

\K N (t) -K hN {t)\ ^0 asN^oo (10) 

uniformly in r = t/N £ [e, (1 — e)e/C\}. 

Our main result is a version of the above theorem that drops the dependence on 
the Conjecture and extends the interval of convergence. 

Theorem 2. For all e > 

\t/N - K hN (t)\ ^0 asN^oo (11) 
uniformly in r = t/N £ [e, e -1 ]. Thus, uniformly in r = t/N £ [e, 1], 

\K N (t) -K hN {t)\ ^0 asN^oo. (12) 
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4. Doubly stochastic matrices and the proof of Theorem [2] 

Definition 3. An entry-wise nonnegative N x N matrix M is called doubly stochastic 
if 

N N 

J2 M i,j = 1 Vj and ^\U;. ; 1 Vz. (13) 

i=l j=l 

The set of all such N x N matrices we denote by DS(iV). DS(iV) is a monoid with 
respect to the matrix multiplication, in particular if A and B are doubly stochastic 
then so is AB. 

It follows directly from the definition that the vector e = (1/N, . . . , 1/iV) is both a 
left and a right eigenvector of any matrix from DS(iV). The corresponding eigenvalue 
1 is the largest (by modulus) one. The following observation shows the importance of 
the asymptotics of the second largest eigenvalue. 

If U is unitary then the matrix M defined by 

M i;j = l^l 2 . (14) 
is doubly stochastic. It is then easy to see that 

N N t 

TrM' = ^(M%= ftl^r. ( 15 ) 

3=1 ii,...,it=l k=l 

and, therefore, the diagonal approximation K^n is equal to 

^W = -^(TrM%, (16) 

where the averaging is performed with respect to the measure induced on doubly 
stochastic matrices by the correspondence (JHJ). In [T2| it was speculated that the 
second largest eigenvalue of such random matrices! M decay, on average, as N- 1 ' 2 . If 
this is indeed so, we immediately get 

N 

(TrM% = 1 + ^(A*) « 1 + (N - l)N~ t/2 -> 1, (17) 

i=2 

for t > 2, which would prove Theorem El While the conjecture of ^2] is still unproven, 
we can use the technique of |T2| to achieve the same result without directly calculating 
the eigenvalues of M. 

Since the eigenvalues of M lie in the unit circle, it is intuitively clear that 
Tr M* = should be a decreasing function in t, at least "on average" . It is therefore 
tempting to show that (Tr M 4 ) n — > 1 for some fixed t as — ■> oo and then conclude 
that this property still holds if we send t to infinity as well. Indeed, this would work if 
the eigenvalues of M were positive real numbers, which is generically not the case for 
matrices M. To circumvent this difficulty, we form a matrix 

A = M T M, (18) 
I an extensive discussion of properties of such random matrices is contained in 
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which is symmetric, doubly stochastic and positive definite. Now we can estimate the 
eigenvalues of M through the eigenvalues of A. 

Proposition 1 ([12]). Second largest eigenvalues o/M and A satisfy 

|A 2 (M)| 2 < A 2 (A). (19) 

Proof. Let v be the eigenvector of M corresponding to A2(M) and orthogonal to e. 
Such eigenvector always exists: if A 2 (M) ^ 1, then we write 

(v, e) = (v, M T e) = (Mv, e) = A(v, e), (20) 

and therefore (v, e) = 0. If A 2 (M) = 1, the geometric multiplicity of 1 is equal to its 
algebraic multiplicity (see, e.g. [Tl]) and we can find v by orthogonalization. 
Then, since all eigenvectors of A are orthogonal, we have 

A 2 (A) = max (Au, u) > (Av, v) = (Mv, Mv) = |A 2 (M)| 2 , (21) 

|u|=l, (u,e)=0 

□ 

Now the idea is to find the power k such that (X(A) k ) N decay sufficiently fast. It 
turns out that k = 3 is enough. 

Lemma 1. The Haar probability measure on the unitary group induces a measure on 
symmetric doubly stochastic matrices A via fli^| ) and hlty . With respect to this measure 
(X 2 (A) 3 ) N = 0(N- 2 ). 

The proof is rather technical and will be given in a later section. This lemma 
enables us to give a simple proof of our main result. 

Proof of Theorem^ We show that (TrM*) — ► 1. If r > e and iV is sufficiently large 
then t = tN > 6, leading to 



|(TrM*) 



N 



,k=2 



N N 



<£|(A,(Mr)|<£<|A,(M)|*) 



k=2 k=2 



< (N - 1) <|A 2 (M)|') < (N - 1) <|A 2 (M)| 6 ) 
<(iV-l)<|A 2 (A)| 3 )=0(iV- 1 ). 

Now we use equation (JTfij) to obtain the sought after conclusion. □ 

5. Proof of Lemma [T] 

We start with the trace of A 3 , 

N 

TrA 3 = l + ^A fe (A) 3 > 1 + A 2 (A) 3 , (22) 

k=2 
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since all eigenvalues of A are non-negative real numbers. Thus we can estimate 
< A 2 (A) 3 < Tr A 3 — 1. Writing out the trace we get 

Tr A 3 = Tr(M T MM T MM T M) 

TV 

= (M T ) lJ (M) J - fc (M T ),, / (M) iifn (M T ) mi „(M) niJ 

i,j,k,l,m,n=l 
N 

i,j,k,l,m,n=l 

Now we average over the unitary group, 

N 

(\{A) 3 ) N + l< £ {\U j ^ j , h U hk U Um U n , m U n ^) N . (23) 

i,j,k,l,m,n=l 

At this point we can apply equation (JBJ with a k = a k and j3 k = b k with 1 < k < 6. The 
indices are a = I, I, n, n) and b = (i, k, k, m, m, i). We denote 

{\U hi U hk U hk U hm U ntm U n ,\ 2 ) = W(t%. (24) 

To cut down the number of the coefficients W that need to be calculated, we list some 
of the properties arising from the symmetries of (|8|). 

(i) Application of a permutation to both a and b leaves W invariant. Due to the 
special form of a and b, the only eligible permutations are shifts: 

= = w kZ (») 

(ii) We can switch around a and b. After that we should apply a right shift to bring 
them back to the special form: a' = (i, i, k, k, m, m) and b' = (n,j,j, I, I, n). Thus, 

Wg£ = K'p- (26) 

Note that the lower indices on the right-hand side are shifted from their original 
order. 

(iii) The coefficients W do not depend on the individual values of j, I and n, but only 
on whether they are pairwise distinct. That is, if 

j = I 4» f = I', I = n<$l' = n' and n = j n' = j', (27) 

then 

W&% = (28) 
The same applies to i, k and m. 

Using these properties we can list all possible values of W. If {j, l,n} are pairwise 
distinct, W/yJJ, = Wl' k 2 '^ and there are N(N — 1)(N — 2) ways to select such j,l,n. 
If I = n but j is different, W-'^ = ^Ikm an d there are N(N — 1) way to choose j 
and n; the cases j = I ^ n and j = n ^ I are analogous. Finally, when j = I = n, 
^i'tm = ^ikm an d there are N ways to choose the value of j. The dependence on 
indices {i, k, m} can be treated in the same way. Table Q lists all contributing terms 
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Contributing terms 


Notation 


Value 


Number of occurrences 


1 O Q 

vv 1,2,3 


^33 


7V°+6iV , '+5iV J — 20iV z — 47V+32 
L>(iV-2)(Ar 2 -l)(7V-l)7V 


N 2 (N — l) 2 (iV — 2) 2 


1 O O O 1 O O O 1 

vv 1,2,3 > rr 1,2,3 ) rr 1,2,3 


^23 


D(N-l)N 


AT2(at- i) 2 (iV-2) 


1 O O 010 OOI 

rr 1,2,2 5 Kr 2,1,2 5 vv 2,2,1 


^22a 


at 7\7"2 i q /vr i a \ 
\>yi\ H-oiV +d ) 

D(N—1)N 


iV 2 (jV- l) 2 


1 O O O 1 O OOI 

vy 2,1,2 ) "2,2,1 5 "1,2,2 


W 22 6 


o\Jy H-oiv+oJ 
D(N-l)N 


iV 2 (iV- l) 2 


" 2,2,1 5 vv 1,2,2 1 "2,1,2 


^22c 


R( N-l-'?'\( M-L~\\ 

D(N-1)N 


AT2(at- l) 2 


rK l,2,3 


Wl3 


8 
D 


iV 2 (iV-l)(iV-2) 


^1,1,1 TT/1,1,1 ^1,1,1 
rk 1,2,2 ' rK 2,l,2 ? vv 2,2,1 


^12 


48 
D 


N 2 (N — 1) 


vv 1,1,1 




720 
D 


iV 2 



Table 1. Factors contributing to (Tr A 3 ). We use notation D = (N + 5)(N + 4) (AT + 
3)(iV + 2)(N + l)N. The remaining terms can be obtained using W/y^ = W^yJ 71 . 
Whenever the first column lists three terms, the last column gives the number of 
occurrences of each of those terms. 



together with the number of times they occur in expansion (|23|). The values of W were 
computed using a program which utilizes formula (|SJ) and the data for Vjv from 
Thus we get 

(Tr A 3 )^ = N 2 (N - l) 2 (iV - 2) 2 w 33 + 2 x 3N 2 (N - 1) 2 (N - 2)w 23 

+3N 2 (N - l) 2 w 2 2a + SN 2 (N - l) 2 w 22b + 3N 2 (N - l) 2 w 22c 
+2 x N 2 (N - l)(N - 2)wj3 + 2 x 3N 2 (N - l)io la + N 2 w u , 

where the factors of 2 arise due to W/'j^ = W^' •'["'. The result is 

/ qvA 3x _ A/ 6 + 16A/ 5 + 105A/ 4 + 370iV 3 + 716A/ 2 + 348A7 - 136 

' 1 (iV + l) 2 (iV + 5){N + 4){N + 3){N + 2) ( ^ 

= l + 5AT 2 + 0(iV™ 3 ), (30) 

which, together with inequality A2(A) 3 < Tr A 3 — 1, finishes the proof. 



6. Discussion 



We presented an alternative and complete proof of the main result of [7j. However 
the main importance of the results of [7j lies in the method of their derivation and in 
establishing a foothold in understand the combinatorics of the form factor expansion. 
Below we formulate two conjectures which can be studied using the methods of jjj- 

The first conjecture is rephrasing an observation made in ^2] regarding the 
asymptotics of the eigenvalues of doubly stochastic matrices M. 

Conjecture 2. For any j > 2 and k (|Aj(M)| fc ) = 0(N~ k / 2 ), where the averaging is 
performed with respect to the measure induced by J7^| ). 
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The second conjecture is a strengthening of our result which would follow if 
Conjecture 121 were proved. It concerns the speed of convergence of the diagonal 
approximation 

Conjecture 3. For any r G (0, oo) 



faster than any power of N . 
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